GLOBAL STABILITY FOR A CLASS OF VIRUS MODELS WITH CTL 
IMMUNE RESPONSE AND ANTIGENIC VARIATION 



MAX O. SOUZA AND JORGE P. ZUBELLI 

Abstract. We study the global stability of a class of models for in-vivo virus dynamics, 
that take into account the CTL immune response and display antigenic variation. This class 
includes a number of models that have been extensively used to model HIV dynamics. We 
show that models in this class are globally asymptotically stable, under mild hypothesis, by 
using appropriate Lyapunov functions. We also characterise the stable equilibrium points for 
the entire biologically relevant parameter range. As a byproduct, we are able to determine 
what is the diversity of the persistent strains. 



1. Introduction 

1.1. Models for in-vivo virus dynamics. A number of population dynamics models have 
been proposed in order to describe the HIV in-vivo dynamics (Perelson & Nelson 1999, Nowak 
& May 2000). Although these models have distinct features, since they attempt to incorporate 
different aspects of the interaction between the virus and the immune system, many of them 
share a common long-term behaviour, evolving towards an isolated equilibrium state (Nowak 
& May 2000). 

The basic model for the HIV in-vivo dynamics is given by a three-by-three, first-order 
system of ordinary differential equations (ODEs) — (Nowak & Bangham 1996, Bonhoeffer, 
et al. 1997, Nowak & May 2000): 



(1) 



x = A — dx — f3xv, 
V = Pxv - ay, 
v = ky — uv. 



In this model, x denotes the uninfected cells, y the infected cells and v the free virus particles. 
The average lifetime of an infected cell is 1/a, while the average lifetime of a virus particle is 
1/u. The total number of virus particles produced by an infected cell is k/a. Healthy cells 
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are infected at a rate f3xv. New CD4+ T cells are produced, in the thymus, at a rate A, and 
die at a rate dx. 

System ([1]) has two equilibrium points: 

(1) The disease free equilibrium: x* = X/d, y* = v* = 0; 

(2) The endemic equilibrium: x* = au/(3k, y* = {(3Xk — dau)/(3ak, v* = {(3Xk — dau)/(3au. 

The long term dynamics of System (pQ) can be entirely described in terms of the dimensionless 
parameter 

0Xk 



(2) 



dau ' 



also known as the basic reproductive ratio. 

If i?o < 1, the disease free equilibrium is a global attractor, and the infection cannot 
persist. If Rq > 1, the endemic equilibrium becomes a global attractor, and the infections 
persists indefinitely. This has been first observed numerically (Nowak & Bangham 1996, 
Bonhoeffer et al. 1997, Nowak & May 2000). Mathematical proofs of these global stability 
characteristics were given by Li & Muldowney (1995), using Hirsch's theory of competitive 
differential systems — see Smith (1995) — and, more recently by de Leenheer & Smith (2003) 
and Korobeinikov (2004a) using a Lyapunov function approach. 

Given the notable ability of the HIV to escape from immune response, there is interest 
in studying models that account for a more detailed immune response as, for instance, the 
role of cytotoxic T lymphocytes (CTLs). An example is the following four-by-four system of 
ODEs (Nowak & Bangham 1996): 



(3) 



x = X — dx — f3xv, 

y = j3xv — ay — pyz, 

v = ky — uv, 

z = cyz — bz. 



System ([3]) extends ([I]) by introducing the z variable, that denotes the CTL response. Infected 
cells are killed at a rate pyz, while antigen stimulation produces CTL cells at a rate cyz. In 
the absence of such a stimulation, CTL cells decay at a rate bz. 

In the same vein, the high mutation rate of HIV naturally leads to the study of the interplay 
between immune response and virus diversity for a number of different strains. The immune 
response produces a selection pressure on these different strains of the virus, as discussed in 
Nowak & Bangham (1996), when studying and numerically analysing a (3n + l)-by-(3n + 1) 
first-order ODE system of the form: 



(4) 



x 
Vi 



X - dx - xYXziPiVi, 

faxvi - aiyi - piyiZi, i = l,...,n, 

k%yi Uii>i, i = 1, . . . , n, 

CiVi^i biZi, i — 1, . . . , Tl. 
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System @ is a slightly generalised form of the system studied by Nowak & Bangham (1996), 
where the more restricted case dj = a, pi = p, ui = u, Ci = c, and bi = b was addressed. 

In all these models, there is the question as whether the long term dynamics approaches 
an equilibrium or, more generally, an attractor, and how this might depend on the initial 
condition. There is compelling numerical evidence — cf. (Nowak &: Bangham 1996, Bonhoeffer 
et al. 1997, Nowak & May 2000) — that Systems §3§ and are globally asymptotically stable. 
However, no mathematical proof of this fact seems to be available. In System there is also 
the question of determining the antigenic diversity at the equilibrium. 

In this work, we study the stability characteristics of the models given by (|4|) following 
a Lyapunov approach. The Lyapunov functional used here has been used before by Kore- 
beinikov & Wake (1999), in the global analysis of three-dimensional predator-prey systems, 
and by Korobeinikov (2004a) and Korobeinikov (2004b) in the global analysis of various virus 
models. More precisely, by using an appropriate linear combination of this Lyapunov func- 
tional, we are able to show global asymptotic stability results for System and hence to 

The plan for this article goes as follows: We close this introductory section with further 
biological background and motivations. In Section [2] we address some preliminary issues such 
as choice of dimensionless variables and parameter reductions. We study the global stability 
characteristics of model ([3]) . In Section O we study the equilibria and global stability of 
model ([!]) under the assumption of unique fitnesses of the strains. In this case, we determine 
the possible equilibria of (|4|), and show than there are 2 n_1 (n + 2) equilibrium points. We also 
show that system is globally asymptotic stable, and we determine what is the global attractor 
in the nonnegative orthant of M 3n+1 . As a byproduct, we characterise the attained diversity 
and show that it is monotonically increasing with the strength of the immune response. Some 
additional results for the case of nonunique fitnesses are also presented. We conclude in 
Section [4] with a discussion of some of the implications of our results. 

1.2. Biological Background and Motivation. The models studied by equation |4] have 
many potential biological applications. Most notably, to within- host infections connected 
to cytotoxic T lymphocytes with antigenic variation including, but not restricted to, HIV 
infection. A better understanding of how the within-host HIV, interacts with immune cells 
seems to be a key factor in the development of effective long-term therapies or possibly pre- 
ventive vaccines for deadly diseases such as the acquired immunodeficiency syndrome (Nowak 
&: May 2000). Mathematical modeling of the underlying biological mechanisms and a good 
understanding of the theoretical implications of such models is crucial in this process. In- 
deed, it helps clarifying and testing assumptions, finding the smallest number of determining 
factors to explain the biological phenomena, and analysing the experimental results (Asquith 
&: Bangham 2003). Furthermore, modeling has already impacted on research at molecu- 
lar level (Nowak & May 2000) and important results have been obtained in modeling the 
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virus dynamics for several infections, such as the HIV (Nowak & Bangham 1996, Perelson, 
et al. 1993, Perelson, et al. 1996), hepatitis B (Marchuk, et al. 1991), hepatitis C (Neumann, 
et al. 1998), and influenza (Bocharov & Romanyukha 1994). 

In the particular case of the HIV infection, the dynamics of the within-host infection goes 
as follows: First, the HIV enters a T cell. Being a retrovirus, once the HIV is inside the T cell, 
it makes a DNA copy of its viral RNA. For this process it requires the reverse transcriptase 
(RT) enzyme. The DNA of the virus is then inserted in the T-cell's DNA. The latter in 
turn will produce viral particles that can bud off the T cell to infect other ones. Before one 
such viral particle leaves the infected cell, it must be equipped with protease, which is an 
enzyme used to cleave a long protein chain. Without protease the virus particle is uncapable 
of infecting other T cells. 

One of the key characteristics of HIV is its extensive genetic variability. In fact, the HIV 
seems to be changing continuously in the course of each infection and typically the virus 
strain that initiates the patient's infection differs from the one found a year ore more after 
the infection. In this respect, the introduction of the different strains in the model is crucial 
for it to be realistic. 

In general terms, one can also say that Model (j4]) is similar in spirit to other models 
such as food-chain models. The latter have attracted substantial interest by a number of 
authors. See for example (Roy & Solimano 1986, Kooi &: Hanegraaf 2001, Kooi, et al. 1998) 
and references therein. However, the presence of more general quadratic terms or of logistic 
terms on the right hand side of the different "strains" leads to a potentially richer dynamics 
than the globally stable present in Model HI See for example (Kooi & Hanegraaf 2001) for a 
bifurcation analysis of certain food chain models. 



2. Preliminaries 

2.1. Parameter reduction. As already noticed in the introduction, a more restricted form 
of (j4|) with a number of parameters being strain-independent has been studied by Nowak 
&: Bangham (1996). It turns out that some parameters in can indeed be taken to be 
strain- independent, as we now show. 

We start by noting that if ki = 0, then V{ decays exponentially with rate U{. Also, if pi = 0, 
then the dynamics of Zi does not impinges on the rest of the system. Thus, without loss of 
generality, we can assume that ki,pi ^ 0, i = 1, . . . , n. In this case, following Pastore (2005), 
we rescale the ViS and fas. In addition, we also rescale the ZjS. More precisely, the change of 
variables 

(5) Vi h-> y-Uj, Zi i-> — Zi and fa — fa 

K Pi Ki 
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takes (SI) into 



(6) 



x = A - dx - x Ya=i Pi v i> 
y i = p iXVi - cnyi - pyiZi, i = l,...,n, 
v\ = kyi-UiVi, i = l,...,n, 

Zi — Ciy%Z{ biZi, % — 1, . . . ? Tl. 

Intuitively, the change of variables ((5|) reflects that only the ratio Pi/ki turns out to be 
important, and that this can be already taken into account in the /3jS, provided we rescale 
the ViS. Moreover, it also shows that the precise value pi does not matter, as long as it is 
nonzero. 

2.2. Dimensionless constants. In Nowak & Bangham (1996), it was already observed that, 
in addition to Rq, the quantity cy*/b is also important in determining the global equilibria. 
In a more precise fashion, Nowak &: May (2000) define 

cdu 

which they term the basic reproductive ratio in the presence of immune response. However, 
we follow Pastore (2005), and find more convenient to write 

where 

ab 

An alternative dimensionless constant is the CTL reproduction number given by 

Po = io(ijW) 
Ro 

Although only two constants among Rq, Iq and Pq are independent, and are sufficient to 
describe the regimes of (0), we have chosen to use both, three constants, as some conditions are 
better characterised by Pq, while much of the algebra in the Lyapunov functional derivatives 
is better handled by expressing them in terms of Rq and 10. Thus, we shall use the strain 
dependent constants: 



2.3. Strain sets. In order to deal with plethora of equilibria that arises in System Q, we 
shall now define some notation for some special set of strain indices. This will allows us to 
deal conveniently with the combinatorial structure of the equilibria. 

Without loss of generality, we shall assume that the strains are indexed by increasing order 

of Rq. 

Let Af = {1, 2, . . . , n}. Then, we define the set of the strong responders as 

S = {i e M\P 1 q > 1}. 
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Definition 1. We shall say that the set S of strong responders is consistent, if 

S = {1, . . . ,m}, 1 < m < n. 

This is certainly the cases, if the Iq satisfy Iq > /q +1 . In particular, this holds if Iq = Iq as 
in the model studied by Nowak & Bangham (1996). 

Given a set of indices T, we define 

Two important definitions are given below: 

Definition 2. We shall say that Z C S is an antigenic set, if 

(8) R > 1 + pi iei 
holds. In addition, if 

(9) R < 1 + pi, it 1- 

also holds, we shall say that X is a stable antigenic set. 

Notice that, if S ^ 0, we have that I = {1} is an antigenic set. Let I be the largest integer 
for which the set J = {1, . . . ,1} is an antigenic set. Then we shall say that J is the maximal 
antigenic set. 

Two important facts about the maximal and stable antigenic sets are collected below: 

Lemma 1. Assume that 5^0, and that the strain basic reproductive numbers are distinct. 

(1) // a stable antigenic set exists, then it is also the maximal antigenic set. In particular, 
stable antigenic sets are unique. 

(2) If M is the maximal antigenic set, then it is a stable antigenic set. 

Proof. (1) Assume that I is a stable antigenic set and let i\ = maxl. If 1 < k < I, then 
ik € T. Indeed, by the increasing ordering and ([8|), we have that 

Ri k >R%>l + P l 

But this contradicts ([9]), thus we must have ik £ T. Now assume that T' = . . . , ii+\} 
is also an antigenic set. Then we must have 

i^ 1 >l + ^ = l + ^ + ^>l + /. 

-'o 

But this again contradicts ([9|) and, hence, that T is stable antigenic. Therefore, it is 
maximal. 

(2) This follows since, in this case, Q cannot be violated. 

□ 
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3. The model with antigenic variation 

In this section we shall study the stability of of ([6]) in the non-negative orthant of ]R 3n+1 
which we shall denote by O. The positive orthant will be denoted by + . 

We observe that the planes z% = and that O are positive invariant sets for ([J]), since the 
field points inwards. 

The equilibria and stability characteristics of ([6]) depend significantly whether the R l s are 
distinct or not. In §3.1 and §3.2, we describe the equilibria and study their stability in the 
case of unique fitnesses, i.e., we assume that if R l = Rq, then i = j. In this case, with the 
adopted order, we have that 

Rl>R^ +1 , i = l,...,n-l. 
Additional remarks when the fitnesses are not unique can be found in Section 3.3. 

3.1. Equilibria. Let J\f = {1,2,... , n}. It turns out that the equilibria of (JH) can be con- 
veniently indexed by (j,J~), where J C M, and either j = or j J. The corresponding 
equilibrium point will be denoted by Xj t j. 
Using this notation, we have 

Lemma 2. System © has 2 n ~ 1 (2 + n) equilibrium points which can be written as 

y „—(^n x ^ nv 1 ^ r>v n ^ n Vl ^ n Vn ai n zi a ' n n Zn 

where 
(1) 

Qo,0 = l' an< ^ Qo$ = Qq,$ = Qo'e = 

(2) If j is such that 1 < j < n, then we have 

QU ( 1% 1 7^7- ^0 = ^0-1 <>"<< Q% o. 



and 



(3) Given J C M , we have 



and 



also 



1 /?* 

i o i o i + Po 
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(4) Given a proper subset J C M , and 1 < f < n,j' J , we have that 



-Kg -Kg -Kg 

for i £ J ', we have 



r>yi - J_ — n Zi — :5lL - 1 

-'O J -Kg 

For i ^ J , and i ^ j' , we have 

Proof. The first equilibrium is trivial. The second type of equilibria is obtained by choosing 
an index j such that Zj = 0, but t)j ^ 0. We can choose only one such j, since this determines 
x. For the other indices i, we set yi = Vi = Zi = 0. The first equation then determines Vj. The 
third type is obtained by choosing a set J of indices, such that, for i € J, we have Z{ ^ 0. 
This readily determines ?/j and Vi. For i J~, we have Hi = Vi = Zi = 0. The first equation, 
then, determines x. Finally, the last equilibria is found by having a set of indices J , as in the 
equilibrium of the third type, and then choosing an index j' J as in the second equilibrium. 
Again, only one such j' can be chosen. □ 

3.2. Stability analysis. We are now ready to study the global stability of the equilibria 
of System [6l Surprisingly, although there is a large number of equilibria, only four of them 
will be globally stable. In what follows, unless otherwise is said, we shall assume that that 
Rq > -R +1 , for £ = 1, ... ,n — 1, and that the set of strong responders is consistent. 

Theorem 1. For system ([6]), defined on O, and with initial condition at its interior, there is 
always a globally asymptotically stable equilibrium given as follows: 

(1) X 0t9 , ifR%<l; 

(2) X lj0J if I < Rl, and Pg 1 < 1. 

(3) If Pg 1 > 1, let J be the maximal antigenic set. Then 

(a) If J is a stable antigenic set, then the equilibrium Xq^j is globally asymptotically 
stable. 

(b) Otherwise, let j' be the smallest integer such that j' J , which exists by virtue 
of Lemma\]\ Then the equilibrium Xy t j is globally asymptotically stable. 

Proof of TheoremUl Following Korobeinikov (2004a), we shall use the following Lyapunov 
function: 

n 

V(x,y,v,z) = x - x*ln(x/x*) + ^2(yi - y*\n(yi/y*)) + 



n n 1 
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where Cj will be a constant to be specified later on. 

Then, using the uniform notation of the the equilibria of Q, that is set in Lemma [21 see 
53.1, we have that 



V =±V(x(t),y(t)Mt)Xt)) 



, n n f)Zi 



i=l 



i= 
it 



P 



dx + 



dx 



(10) 



a -of j— -dkf^ + E w N - * - 



+ 



1 



O yi - — 

^hJ Ti 



i=l i=l i=l 

For the first two equilibria, we consider the Lyapunov function (|10p . with Cj = a^jk. Then, 
on using the structure of equilibria of ([6]), we may write V as follows: 



V =A 



n n s^Vi n p>Zi 



i=l 



a,; 



=1 J o 



A 2 QJ 



dx 



i=l 



i=l 



i=l 



■in 



For X Oj 0, we find, using Lemma [21 that 



V = 2A 



f/.C 



Since -Rq < 1, for i = 1, . . . , n, and 



i=l 



A 2 

dx + — > 2A, 
dx 



i=i atI o 



we have that V < in + . Hence, that X is globally asymptotically stable in this case. 
Now, suppose that 1 < Rq, and that Pq < 1. In this case, Lemma [2] yields that 



F =A 



X 2 



dx + — 1 



R^dx 



01 



1 



1 \ xt>i aid 1 



- 1)^+ 



i=l 



+ pA 



ai 



1 



1 1 

^0 J o 



n 

^ 5-/ n . Ti ' 

i=2 ° li 



I'l 



The last term in V is clearly negative. We also observe that, since R l < Rq, for 1 < i < n, 
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Also, since Pq 1 < 1, then we have that Rq < 1 + Rq/Iq. Therefore, the last three terms in the 
expression for V are negative. 

For the remaining terms, let us write 



Rl 



X 



+ 1 



1 



X 



Then we have that 



dx + 



R^dx 



> 2 



and that 



R\dx 



Rl) a x V Rl) yi /3i v > v l ~ V Rl) 



Thus V < in + , and hence we have that Xi g is a globally asymptotically stable equilibrium. 

Finally, if Pq > 1, then let J be the maximal antigenic set. First, we assume that J is 
stable antigenic and show that Xq^j is globally asymptotically stable. In this case, we use 
the Lyapunov function (fT0|) . with Cj = x*(5i/ui. 

Using Lemma El this can be further recast as V = V\ + V2, where 



Vi =A 



1 




3 T 

1 + Po. 





dx + 



A 



dx (l + p|f) 



>E 



/3i xvi 



Xk 



E 



R Vi. 



V2=^2 aiVi 



R'n 



1 



t—i n ■ T l 



.1 + Po 7 

We treat V2 first. The last term is clearly negative. Also, since for i J, we have that 

R'n 



1 + Po 7 



< 1 and thus that yj 



*0J 



i? 



1 + Po 7 



1 



< 0. 



Therefore, V2 < 0, when J 7^ A/". 
Let 



Po 7 



1 + Po 7 



and rji 



1 + 



% G J. 



Then, we may write Vi as 





1 




V\ = X 


3 - 





A 



E 

ie.7 



A 2 ^ 



^(1 + Po) 
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For each i 6 J, we have 



AS 



and that 



A 2 



cfe (1 + 



^ 2 



> 2A- 



1 + 



After combining these estimates and summing for i £ J, we get that V\ < and thus we 
have the result, If J7" is a proper subset of M. In the case that J = A/", we have V < 0, with 
equality ocurring only when 



A 



x = —On q and — = — . 



Inasmuch this plane is not invariant by the corresponding flow — other than the point Xq j — 
we have global stability as a consequence of LaSalle's theorem (LaSalle 1964). For the fourth 
point, we use a mix of the two Lyapunov functions above, namely: 



V(x,y,v,z) =x - x*ln(x/x*) + ^ (y* - y*ln(yi/y*)) + 

i=l 

n „ n 

+ x *J2~(vi- v* ln(VO) +PE - & ~ < M^/zl) + 

— ^ it - £ — J r- 



i=l 



+ ~t( v f -Vj>Hvj>/v*,)) . 
Computing V and using the uniform notation, we find that 



V =A 



i + Q f,j + E dp, j + q}.j E Qf j + E% £ + " 

i=i J o 



i=l 



i=l 



^0 



x2 Qj',J 
dx 



A E 



Pi^ji^j XVi , , „^ v V V 



i=l 



ai yi 



n f^Vi n 

i=l % 1 1=1 



+ 



i=l 



O yi - — 
J 



On using Lemma [2 and that 
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where each term in the sum is positive, but smaller than one, we rewrite it as V = V\ + V2 + 
V3 + V4, where 



Vo 



1 



2X_ 



1 



R R Q 



dx + 



3A 



A 2 \(3ji xv ji daj'R J Q yy 
dxR% af Vj> Pj> v f 



A 2 



dx{R%) 2 



2 n J 



AVo J 



pL_ 

R% dx(R^) 2 



A y^ Pi xvj _ \k_ Rq yj 



Va 



pXzj 



1 



„v „-/ 
n -"-0 1 



The terms Vi, V2 and V3 can be treated similarly as in the previous equilibria and are all 
nonpositive in the interior. 

As for V4, first we observe that 



1 



if j' S, then we that 



oi' S' r i' ~ oi' ( ^ 1 ^ r° 



1 



1 



R'n 



T)J T)J T. 



^o'-l-Po r -^<<-l-^<0 

in Jn 



If j' € 5, then let J' = J U {j'}. Then we have that 

l_J__tfL_J_-J_f jH 



T>J T>J TJ 



l 



Since J is the maximal antigenic set, we must have that 

r{ - 1 - Pcf < 0. 



□ 



3.3. Nonunique Fitness. Given the non-generic nature of this case, we shall only briefly 
discuss the stability when some of the strains have the same fitness, i.e., there exists at lest 
one index set T, such that R l = R J , for i, j 6 T. Notice that, in this case, we have non-isolated 
equilibria. 

We start by observing that the computation with the Lyapunov function for the equilibrium 
Xq does not depend on the uniqueness of fitness. Hence we have 

Corollary 1. If R l < 1, for i = 1, . . . ,n, then Xq$ is a globally asymptotically stable equi- 
librium. 



The case ij = Iq, 1 < Rq and Pq 1 < 1 can also be partially treated: 
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Proposition 1. Let V be the set of indices % G M such that R = Rq. Let E^ be the set 

satisfying 

A . , „ \ „ dx* — A 



x 



dR% 



Hi 



Vi 



o, j#r, 



X" 



Vj > 0, 



k 



u 



y h and z { 



eAf. 



7/x(t,xo) is a solution of with initial condition xq, then 

x(i) — > Er, as t 



oo. 



Proof. Let us denote the omega set of x(£,xq) by f2(xo). As shown in (Pastore 2005), the 
solutions to ([!]) are bounded in R+ n+1 . Hence, ri(x ) is compact. 

Using the same Lyapunov function for X 1 q as in Section 3.2, we find that 

A 2 



V =A 



Rq 



dx + 



Rq- 



R^dx 



A 



Pi 1 



XV\ 



,yi 



RoJ vi 



Pi Ul 



1 

R 



+ -£• 



i=l 



1 1 







The same calculations in section 3.2 shows that V < 0. However, notice that V = in E-p. If 
s € 1^ and S C 1^ is closed, then let 

d(s, S) = min lis — s'll. 

s'£S 

LaSalle's invariant principle then yields that, 



lim d(x(t),E T ) = 0. 

t— >oo 



□ 



Corollary 2. If Rq > i?Q for i = 2, . . . , n, then Xiq is a globally asymptotically stable 
equilibrium, when 1 < Rq < 1 + Rq/Iq. 

We have performed numerical calculations of ([!]), using a high order Runge-Kutta method, 
which suggest that in the case treated by proposition 1, a solution of System @ will converge 
to a unique equilibrium point in Er, that depends only on the initial condition. 

Finally, when the viable set of strains is not the full antigenic variation, we have 

Proposition 2. If Iq = Iq, Pq > 1. Assume that J ^ X is a stable antigenic set, then we 
have that Xqj is a globally asymptotic stable equilibrium. 

Proof. In this case, the estimate V\ < remains valid. Moreover, if J ^ X, then we must 
have V2 < 0, which yields the result. □ 
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4. Conclusions 

In this work, we have performed a thorough study of System (JU). As a preliminary result, 
we have shown that, when the both the virus production rate and the CTL interaction 
rate are nonzero for all strains, then System ([3]) is dynamically equivalent to System ([6]), 
which has strain-independent virus production and CTL interaction rates. In particular, the 
precise nonzero values of the CTL interaction rates are completely irrelevant to the dynamical 
behaviour of the system. This seems to suggest that a more refined model which is able to 
capture this difference is needed. 

We have also identified all the 2 n_1 (2 + n) equilibria of (J6j) in Lemma [2j When n is 
large this can be quite a large number. Nevertheless, under the hypothesis of unique fitness, 
we were able to show that only four of them are dynamically relevant. More precisely, we 
assume that the virus basic reproduction rates, Rq are distinct and that the CTL reproduction 
numbers, Pq have the same ordering as the R . This last condition is automatically satisfied 
if Iq = Iq for all i. In this case, Theorem [1] shows that, if the largest reproduction number, 
Rq, is smaller than one, the the disease- free equilibrium — X g in the notation of Lemma[2] — is 
globally asymptotic stable. In this case, no strain is viable and the infection dies out. On 
the other hand, ii Rq > 1, but P < 1, then only the first strain survives, and the infection 
persists. If Pq > 1, then we have that the two outcomes are possible: either a (unique) 
stable antigenic set exists and then Xq j is globally asymptotic stable. In this case, the set 
J determines the antigenic diversity. In other words, a strong immune response generates a 
larger antigenic variation. Alternatively, there exits a pair (j' , J) such that the point Xy j is 
globally asymptotic stable. In this case, the strain with the weakest fitness will not actually 
trigger the CTL response at all in the long run. In the case of absence of antigenic variation, 
i.e. System ([3]), then only the first outcome is possible. We were unable to interpret in a 
biological sense the combinatorial conditions of existence of a stable antigenic set, and we 
believe that this should be addressed in the future. The results presented in Theorem [1] show 
rigorously some of the inferences that have already been made in Nowak & Bangham (1996) 
based on extensive simulations of System ([5]). 

We have also shown some results for very special cases in which the R s are not distinct. 
In these cases, the equilibria is not isolated and this complicates the matters further. Also, 
we have not addressed that case when the set of strong responders is not consistent, and this 
might also merit further study in the future. 
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